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Abstract 

m 

This paper deals with the dynamical system that generalizes the MIC-Kepler 
system. It is shown that the Schrodinger equation for this generalized MIC-Kepler 
system can be separated in spherical and parabolic coordinates. The spectral prob- 
lem in spherical and parabolic coordinates is solved. 
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1 Introduction 
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The system described by the Hamiltonian 
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where c\ and C2 nonnegative constants, later on will be called the generalized MIC-Kepler 
system. (We use the system of units for which h = m = e = c= 1.) 

The MIC-Kepler integrable system was constructed by Zwanziger pQ and rediscovered 
by Mcintosh and Cisneros [2|. This system is described by the Hamiltonian 



^o = ^HV- S A) 2 + ^-i (1.2) 



where 



1 r 

A = — -(y, — x, 0), and rotA = — . 

r(r — z) r 6 

Its distinctive peculiarity is the Coulomb hidden symmetry given by the following 
constants of motion (jl.2j) 



(-iV - sA) x J — J x (-iV - sA)] + ^, J = r x (-iV - sA) - s-. (1.3) 



Here, the operator J defines the angular momentum of the system, while operator I 
is the analog of the Runge-Lenz vector. These constants of motion, together with the 
Hamiltonian, form the quadratic symmetry algebra of the Coulomb problem. For fixed 
negative energy values the motion integrals make up algebra so(4), whereas for positive 
energy values - so(3.1). Due to the hidden symmetry the MIC-Kepler problem is factorized 
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not only in the spherical but parabolic coordinates as well. Hence, the MIC-Kepler system 
is a natural generalization of the Coulomb problem in the presence of Dirac's monopole. 
In both cases the monopole number s satisfies the Dirac's rule of charge quantization 
s = 0,±1/2,±1,.... 

The MIC-Kepler system could be constructed by the reduction of the four- dimensional 
isotropic oscillator by the use of the so-called Kustaanheimo-Stiefel transformation both 
on classical and quantum mechanical levels |3] . In the similar way, reducing the two- and 
eight- dimensional isotropic oscillator, one can obtain the two- |3] and five- dimensional 
[Hj analogs of MIC-Kepler system. An infinitely thin solenoid providing the system by 
the spin 1/2, plays the role of monopole in two-dimensional case, whereas in the five- 
dimensional case this role is performed by the 577(2) Yang monopole endowing the 
system by the isospin. All the above-mentioned systems have Coulomb symmetries and 
are solved in spherical and parabolic coordinates both in discrete and continuous parts of 
energy spectra [7] . There are generalizations of MIC-Kepler systems on three-dimensional 
sphere jH| and hyperboloid [Oj as well. The MIC-Kepler system has been worked out from 
different points of view in Refs. [HI HU E21 ECU HU . 

For integer values s the MIC-Kepler system describes the relative motion of the two 
Dirac's dyons (charged magnetic monopoles), where vector r determines the position of 
the second dyon with respect to the first one £Q. For half-integer s the presence of the 
solenoid magnetic field, endowing the system with the spin 1/2, is presupposed (see, e.g. 

EU) 

The Hamiltonian (jl.lj) for s = and q^O (i — 1, 2) reduces to the Hamiltonian 

* = -| A - i + T i ^ + T i S' ' L4 > 

2 r r[r + z) r[r — z) 

of the generalized Kepler-Coulomb system [TH] . 
The potential 

V= + —, r + r, (1.5) 

r r[r + z) r[r — z) 

is one of the Smorodinsky-Winternitz type potentials [16j. The Smorodinsky-Winternitz 
type potentials where revived and investigated in the 1990 by Evans ^7j. In the case 
where C\ = C2, the potential (jl.5|) reduces to the Hartmann potential that has been used 
for describing axially symmetric systems like ring-shaped molecules |18 j and investigated 
from different points of view in Refs. JSj-jSI]- In particular, the (quantum mechanical) 
discrete spectrum for the for the generalized Kepler- Coulomb system (|1.4jl is well known 
j2H 123 l2~9*j . even for the so-called (q,p) -analogue of this system [21]. Furthermore, a 
path integral treatment of the potential (|1.5[) has been given in Refs. [2SII2ZI- Recently, 
the dynamical symmetry of the generalized Kepler-Coulomb system has been studied in 
Refs. [2Hl EDI EI], ^ ne classical motion of a particle moving in the potential has 
been considered in jHHj, and the coefficients connecting the parabolic and spherical bases 
have been identified in [21, as Clebsch-Gordan coefficients of the pseudo-unitary group 
577(1,1). 

The purpose of the present paper is to further study the bound states of the generalized 
MIC-Kepler system in spherical and parabolic coordinates. 
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2 Spherical Basis 

The Schrodinger equation with Hamiltonian (jl.lj) in spherical coordinates (r, 9, tp) may 
be solved by seeking a wavefunction ip of the form 

if>{r,e,<p)=R(r)Z(6,<p). (2.1) 

This amounts to finding the eigenfunctions of the set {H, J z , M} of commuting operators, 
where the constant of motion M reads 

1 + cos 9 1 — cos 9 

Here J 2 is the square of the angular momentum, J z = s — id/ dip its ^-component and 
J z ip = mip. 

After substitution the expression (j2.1j) the variables in the Schrodinger equation are 
separated and we arrive at the following system of coupled differential equations: 

1 d ( . 8Z\ 1 ( d 2 \„ 

sm9— + — - 4c x \Z + 



sin9d9 \ 89 J 4cos 2 f \d<p 2 



4 sin 2 | 



-AZ, (2.3) 



r z ar \ ar J r A \ r J 

where A is a separation constant in spherical coordinates. 
The solution of (j2.3J) is easily found to be 

/ fj\ mi / f)\ m2 

zM(e,w5 1 ,5 2 ) = N jm (5 1 ,62) cos- fsin|j P^\<x*0)e**->* , (2.5) 



where mi = \m — s\ + S\ = y(m — s) 2 + 4ci, m 2 = \m + s\ + 82 = y (m + s) 2 + 4c2, 
m + = (|m + s\ + |m — s|)/2 and P^ a ' 6 ) denotes a Jacobi polynomial. The quantum 
numbers m and j run through values: m = —j, + — and 

\m + s\ + \m — s\ \m + s\ + \m — s\ 
J = g , g + 1,.... 

The quantum numbers j, m characterize the total momentum of the system and its 
projection on the axis z. For the (half)integer s j,m are (half) integers. 
Furthermore, the separation constant A is quantized as 

( . 5 X + 5 2 \ ( . 5 X + 5 2 A (n nX 
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The normalization constant Nj m (5i, 8 2 ) in (J2.5)) is given (up to a phase factor) by 



Nj m (5i, 8 2 



(2j + 6 1 + 8 2 + - m+)\T{j + m + + 8 t + 8 2 + 1) 

2 m+ ^ 2 s ^+ 2 nT(j - m_ + 5i + l)T(j + m_ + 5 2 + 1) ' 



(2.7) 



where m_ = (|m + s| — \m — s\)/2. The angular wavefunctions Zj„ [see Eq. (J2.5|) ] are con- 
venient to call the ring-shaped monopole harmonics by analogy with the term "monopole 
harmonics" studied by Tamm |32]- These ring-shaped monopole harmonics generalize the 
functions studied by Hartmann [TH] in the case s = 0, 8\ = 8 2 . Due to the connecting 
formula 



A+2) C^(x) = (2X) n P, 



(x 



(21 



between the Jacobi polynomial Pjf 1 '^ and the Gegenbauer polynomial C*, the case s = 0, 
8\ = 8 2 = 8 yields 



^,^^) = 2M + »r( H+ ^i)J^ 



(2j + 28 + l)(j- 


m\)\ 


4vr 2 r(j + 


m 


+ 25+1) 



smi 



\rn\+S c \rn\+6+^, d)em ,- 



3-\m\ 



(2.9) 



the result already obtained in Ref. |15j . [In (J2.8)) (a) n stands for a Pochhammer symbol.] 
The case 8 = (i.e. c\ = c 2 = 0) can be treated by using the connecting formula 



PPix) 



(-2)1 



7T 









\m\ - 




(1 - **) 









J j-\m\ 



(2.10) 



between the Gegenbauer polynomial and the associated Legandre function (23]. In 
fact for the 8 = 0, Eq. (j2.9j) can be reduced to 



zS(0,¥>;O,o) 



^' + 1 )O'-H)^H (cos0)e ^ ; 



\ 4vr(j + \m\)\ 
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an expression (up to a phase factor) that coincides with the usual (surface) spherical 
harmonics Yi m (0,ip). 

Let us go now to radial equation ()2.4j) . The introduction of ()2.6|) into the (|2.4|1 leads 

to 



Id/ 2 dR 



1 



r 2 dr \ (ir / r 2 \ 



8 X + 8 2 



3 + 



8i + 8 2 



+ UR + 2(E+~ > )R = 0, (2.12) 



which is reminiscent of the radial equation for the hydrogen atom except that the orbital 
quantum number I is replaced here by j + (81 + 8 2 )/2. The solution of ()2.12j) for the 
discrete spectrum is 

Rn](r) = Crvfa, 8 2 )(2ery +h ^e- £r F (-n +j + 1; 2j + ^ + 5 2 + 2; 2er) , (2.13) 



where n=|s| + l,|s| + 2,....In (|2.13|) . the normalization factor C n j(8i, 82) reads 

2e 2 



C n j(8i, 8 2 ) 



r (n + j + 81 + 8 2 + 1) 



T(2j + 8 1 + 8 2 + l)\ (n-j-1)! 
and the parameter e is defined by 



n + 



61+62 



The eigenvalues E are then given by 

E = 



2 (n+ ^ 



2 ' 



(2.14) 



(2.15) 



(2.16) 



In the limiting case 8\ = 82 = 0, we recover the familiar results for charge-dyon bound 
system JTj. 



3 Parabolic Basis 

Let us consider the generalized MIC-Kepler system in the parabolic coordinates. In the 
parabolic coordinates £, r\ e [0, 00), if G [0, 2-71"), defined by the formulae 



(3.1) 



x = y£r]cos(p, y = ^£77 sin<£>, z = -(£-r)), 



the differential elements of length and volume read 



while the Laplace operator looks like 
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(3.2) 



A 



d_ t d_\ + d_(d_ 
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1 d 2 



if] dp 
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The substitution 



i(m—s)ip 



2ll 



(3.3) 



(3.4) 



separates the variables in the Schrodinger equation and we arrive at the following system 
of equations 



^ ) 

d ( rf$2 X 



drj \ drj 



+ 



E m\ 1 1 



E ml l n I 

—n - -13+ - 

2 ' Ar] 2 H 2 



$1 = 0, 



$5 



(3.5) 
(3.6) 



where /3 — is the separation constant. 

These equations are analogous with the equations of the hydrogen atom in the parabolic 
coordinates [33]. Thus, we get 



i(m—s)ip 



271 



where 



Here n\ and n 2 are nonnegative integers 
\m - s\ +Si + 1 3+1 

Tli = 1 , 

2 2e 



T(ni + mi + l) _i2 s 

— e 2 (ex) 2 F(— rii) ra; + 1; ex) 



n 2 



\m + s\+5 2 + l (3-1 



2s 



(3.7) 



(3.8 



(3.9) 



From the last relations, taking into account (j2.16|) . we get that the parabolic quantum 
numbers n\ and n 2 are connected with the principal quantum number n as follows 



n = n 1 +n 2 + 



\m — s\ + \m + s\ 



+ 1. 



(3.10) 



Excluding the energy E from Eqs. ()3.5|) and (|3.6|) . we obtain the additional integral of 
motion 



X 



did 



d Id 



2 €-r) 2cm 



£ — n d 2 . f 2 + ii 2 d 

+ is 

2£r] dip 2 £v(£ + v)dy> 



(3.11) 



— s 



+ 



with the eigenvalues 

3 — e I ni — n 2 + 



in — s 




m + s 


+ S X -S 2 \ 


2 1 



(3.12) 



and eigenfunctions ^ { n % 2m (£, ij, (p; Si, S 2 ). 

In Cartesian coordinates, the operator X can be rewritten as 



X 



d 2 d 2 



d 2 



d 2 



X 



+ is- 



r + z 



d d 
x- — y- 



dx 2 dy 2 J dxdz dydz r(r — z) \ dy dx 



d 2 r + z r — z r + z z 

~74 s ~i \ + Cl ~t — I — V ~~ ° 2 ~t V + ~> 

dz r[r — z) r[r + z) r[r — z) r 



(3.13) 



so that it immediately follows that X is connected to the z-component f z of the analog 
of the Runge-Lenz vector (jl.3j) via 



- f r — z r + z 

X = I z + Cl— ; r - C 2 - 



r(r + z) r(r — z) 



(3.14) 
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and coincides with I z when c\ = Ci = 0. 

Thus we have solved the spectral problem in spherical 



H^ = E^, M^= (j + ± Y ± j U + ^y^ + 1 J^' -U = m^ (3.15) 

and in parabolic coordinates 

Hip = Eip, Xip = pip, j s i() = mip, (3.16) 

where H, J z , M and X are defined by the expressions ([□}, (ESI), and (l3~T4l . 

It is mentioned that all the formulae obtained for s = yield the corresponding 
formulae for the generalized Kepler- Coulomb system [TH] . 
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